The correlation exponent Kp of the one-dimensional Kondo lattice model 
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We present results for the correlation exponent Kp of the Tomonaga-Luttinger liquid description 
of the one-dimensional Kondo lattice as a function of conduction electron density and coupling 
constant. It is obtained from the first derivative of the Fourier transform of the charge-charge 
correlation function. We also show that the spin correlation function can only be described in this 
picture if we include logarithmic corrections, a feature that had been previously overlooked. A 
consistent description of both charge and spin sectors is then obtained. Finally, we show evidence 
that the spin sector of the dimerized phase at quarter-filling is gapless. 
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The Kondo lattice model is the simplest model be- 
lieved to describe the low energy physics of heavy fermion 
materials ii Its one-dimensional version has been thor- 
oughly studied in the last 10 years and a great deal 
of understanding has been gained. However, some out- 
standing issues remain, some of which may have impli- 
cations in the higher dimensional cases. For example, 
the question of whether the localized spins should be 
counted in a Luttinger's theorem determination of the 
size of the Fermi momentum is still controversial i^'^'^'^ 
Furthermore, even the phase diagram still presents some 
novel surprising phases: at quarter conduction electron 
fining the spins are dimerized and the charge sector is 
gapped.'' The latter phase may be at the origin of the 
spin-Peierls phase observed in the quasi-one-dimensional 
organic compounds (Per)2M(mnt)2 (M=Pt, Pd)i^ At a 
generic incommensurate filling however, the system is 
gapless in both the spin and charge sectors and it is 
reasonable to assume^ that it is a Tomonaga-Luttinger 
liquid (TLL)iifi In this paper, we will assume that this 
is the case. In an attempt to systematically character- 
ize this behavior, we have determined the non-universal 
TLL exponent Kp as a function of coupling constant and 
conduction electron density. We found that a consistent 
picture of charge and spin sectors can be obtained, only 
if logarithmic corrections are included in the spin corre- 
lations. Moreover, we also show that the spin excitation 
spectrum of the quarter-filled case is gapless. We give ar- 
guments showing that the presence of dimerization and 
the absence of a spin gap are not mutually exclusive. 

We considered the one-dimensional spin-^ Kondo lat- 
tice Hamiltonian with L sites 
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where Cjcr annihilates a conduction electron in site j with 
spin projection (t/2, is a localized spin-i operator and 



2 So 



^cTapc- g is the conduction electron spin 



density operator. J > is the Kondo coupling constant 
between the conduction electrons and the local moments 
and the hopping amplitude has been set to unity to fix 



the energy scale. We studied the model with the density 
matrix renormalization group (DMRG) techniqueii with 
open boundary conditions. We used the finite-size algo- 
rithm for sizes up to L = 120 keeping up to m = 600 
states per block. The discarded weight was typically 
about 10~^ — 10~^ in the final sweep. 

TLL's with periodic boundary conditions and SU(2) 
symmetry have charge and spin correlation functions 
given asymptotically hy^ 
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where 6n (x) = n{x) — {n{x)), S'^{j) = Sj +Sj, Kp is the 
non-universal charge correlation exponent and kp is the 
Fermi momentum. Local charge perturbations, such as 
introduced by impurities or boundaries, induce density 
oscillations, called Friedel oscillations. In the case of a 
TLL, they take the formiiiiS 



{6n (x)) = a 
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The main goal of this work is to present Kp as a func- 
tion of the conduction electron density and the Kondo 
coupHng J for the one-dimensional Kondo lattice model. 
A previous work^ determined Kp, but only for the den- 
sity n = 2/3. Besides, in that work, the authors argued 
that the system has a "large" Fermi surface, with 2kp = 
7r(n + 1) (mod 2tt), not a "small" one with 2kp — nn 
(mod 2tt). Indeed, under some assumptions, the pres- 
ence of low-lying excitations with momentum 2kp can be 
provedj^i Assuming a "large" Fermi surface, the numeri- 
cal results show that the dominant term in the charge 
Friedel oscillations is the second one in Eq. (EJ. From 
the decay of the envelope function of this term, Kp was 
determined at n = 2/3!^ However, more recent work has 
called into question the presence of a "large" Fermi sur- 
face, particularly for small If the Fermi surface is 
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Figure 1: Charge gap vs l/L for J = 10 and densities n — 
0.4 and n = 0.6. Here, E{N) is the ground state energy 
of the sector with A'' electrons. The dashed lines are fits to 



J=0.35 
J=0.5 




Figure 2: Fourier transform C(g) versus momentum for sev- 
erals values of J, L = 40, and density n — 0.8. The arrows 
indicate the position of the cusp. 



small, both terms in Eq. (pj oscillate with the same pe- 
riod at n = 2/3 and the envelope function method cannot 
be unambiguously applied. In order to avoid this ambi- 
guity, we determined Kp from the first term in Eq. 1^, 
or equivalently, from the derivative of the Fourier Trans- 
form of the charge-charge correlation function at g = 



Kp = Tr 



dC{q) 



(4) 
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where 



j.k 

This method has been shown to give very accurate re- 
sults by Daul and Noackii^ These authors determined 
the exponent Kp (by the DMRG technique) for the one- 
dimensional Hubbard model and found good agreement 
with the exact results. For this reason, in the present 
work we will use this procedure to estimate the exponent 
Kp. 

Haldane has conjectured that the TLL is the generic 
universality class of one-dimensional gapless systems ~i 
Although a rigorous proof usually relies on the integrabil- 
ity of the model, renormalization group arguments con- 
firm this conjecture in paramagnetic phasesiitt Much less 
is known about the case of systems with ferromagnetic 
ground statesi^^ However, even in this case, the spin sec- 
tor usually decouples from the charge sector and it is pos- 
sible for the latter to remain a TLL. The one-dimensional 
Kondo lattice model is ferromagnetic for sufficiently large 
Jpi& In Fig. n] we show its charge gap as a function of sys- 
tem size for J = 10 (inside the ferromagnetic phase) and 
the densities n — 0.4 and n = 0.6. The extrapolated val- 
ues suggest that the ferromagnetic phase, Hke the param- 
agnetic one, has no charge gap. Thus, it is quite natural 



to expect that, inside the ferromagnetic phase, the charge 
sector may also be described as a TLL and we will assume 
this to be the case. In fact, as we will see below, our re- 
sults are consistent with this assumption. Note that this 
appears to happen also in the case of the Hubbard model 
with next-nearest neighbor hoppingiii 

We first focus on the general g-dependence of C (q) . In 
Fig. 121 we present the Fourier transform of the charge- 
charge correlation function for n = 0.8, L — 40, and 
several values of J. We have checked that the qualitative 
behavior of C{q) presents no finite size effects, and also 
observed that the simple sum rule C(0) = is satisfied 
(within the accuracy of the DMRG) for all values of den- 
sity and Kondo coupling J shown. For small values of 
J and all densities, C{q) increases linearly with q up to 
q = nn, and then saturates at C(g) = n for n < q/ir < 1. 
On the other hand, for large Kondo coupling C{q) in- 
creases linearly with g up to g = 2?™ (mod 27r) < tt and 
then saturates at C(g) = n (l — n) for n < 1/2 (n > 1/2) 
and 27rn (mod 27r) < g < tt. 

In order to get some insight into the behavior of C(g) 
we consider free fermions with spin-S in a one dimen- 
sional nearest-neighbor tight-binding lattice. In this case, 
the Fourier Transform of the charge-charge correlation 
function Cq {q) is 



C^(g) _ r g/27r < g/vr < 2m, 
(25" +1)^1 m 2m < q/n < 1, 



(5) 



where m = min [n/(2S' + 1), 1 - 71/(25" + 1)]. We will 
need two particular cases, with the restriction ti < 1. 
For spin-i fermions 



g/TT < g/TT < n, 
n n < g/vr < 1, 
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Figure 3: The exponent Kp obtained with Eq. Elas a function 
of J with L — 40. The densities are indicated. Inset: Kp vs 
density for J=0.35 and J=0.5. 

while for spinless fermions, if n < 1/2, 

r^ln\ - / ^/^^ < g/TT < 2n, 
n 2n<qlTT<l, 

and if n > 1/2, 

ror^^-j '^/^^ < g/^ < 2(1 - n), 
^o^'^^~\l-n 2(1 - 7i) < g/vr < 1. 

Our results for C{q) in the one-dimensional Kondo lat- 
tice model all tend to the free spin-^ case when J ^ 
(see Fig. El, as expected. Besides, for J ^ 0, C(g) tends 
to the 5 = case, Cq(ci). This is also to be expected, 
since in this case the conduction electrons form unbreak- 
able mobile singlets with the locaHzed spins, effectively 
behaving like spinless fermions fl^ 

The cusp of C{q) &i q — 2k p = mr, for small values 
of J is the signature that the charge density oscillations 
(Eq.^ are dominated by the 2kp term. As we increase J 
the cusp moves to g = 27rn (mod 2tt). At first sight, this 
might seem like an indication that the system crosses over 
from a 2fc|l-dominated region to a 4fc|^-dominated one as 
J increases. Indeed, this is what happens in the Hub- 
bard model when we increase the on-site repulsion UM- 
However, since Akf, — 27rn = Akp (mod 2tt), the change 
might be due to a phase transition from a small Fermi 
surface to a large onei^ Indeed, there have been indica- 
tions of an intervening ferromagnetic phase at interme- 
diate values of J^^ which could give rise to this change. 
Unfortunately, the study of other (spin-spin) correlation 
functions has not shed any light on the issues^- the size 
of the Fermi surface for intermediate values of J remains 
an open question. 

In Fig. we show the exponent Kp calculated through 
Eq. Q for several values of n and Kondo coupHng J. For 
all densities we see that Kp tends to unity when J ^ 0, 



Figure 4: The exponent Kp vs 1/L ioi J = 0.35, n = 0.4 
(circles) and n — 0.8 (squares). The dashed lines are fits 
according to Eq. H with Kp = 0.55 (0.76) a = 17.0 (10.3) 
and b = -196.1 (-120.5) for n = 0.4 (0.8). 



in agreement with our expectation that the system tends 
to a non-interacting spin-^ electron gas with Kp — 1. On 
the other hand, in the strong coupling limit Kp ~ 0.5, as 
expected for free spinless fermions (cf. Eqs.01and[2||. We 
have observed that for densities n > 1/2 and J ~ 1.5, Kp 
attains its smallest values. This is a region where charge 
oscillations are enhanced (see, e. g.. Fig. 1 of Ref.|3). Ac- 
tually, for some values of J and n > 1/2 we were unable 
to determine Kp, as for example at J = 1.7 and n = 0.9. 
At these points, our data did not satisfy the simple sum 
rule C(0) = even increasing the truncation m up to 
m = 600. The exponent Kp also did not converge as a 
function of to. A truncation of to = 600 is more than 
enough to get precise values in other parameter regions. 
Typically, for small values of J and L — AO, Kp obtained 
with truncations m = 400 and m — 600 differ by less than 
5.10-2 and C(0) ~ 10"* with to = 400. It is interest- 
ing to note that this region where the charge oscillations 
are strongest corresponds to the ferromagnetic phase at 
intermediate JM For completeness, we also show in the 
inset of Fig. El the dependence of Kp on the density for 
J = 0.35 and J = 0.5 with L = 40. As we can see, Kp 
decreases with increasing J. 

Given this qualitative behavior of Kp as a function of 
J, we now set out to determine some quantitative values 
of the exponent. For this, we must be careful to take 
into account finite-size effects. In Fig. ^ we show Kp for 
the densities n = 0.4 (circles) and n = 0.8 (squares) as a 
function of 1/L at J = 0.35. In order to incorporate the 
finite-size dependence, we determined the extrapolated 
exponent assuming that Kp{L) behaves Hke 

Kp{L) ^ Kp + a/L + b/L^. (9) 

In Fig.^the dashed fines are fits to our data using Eq. El 
The exponents Kp obtained through the fits are 0.55 and 
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Figure 5: The spin-spin correlation function for densities 
n = 0.4 (a) and n = 0.8 (b) (J = 0.35 and L = 120). The 
dashed line is a fit of Eq. ED with a — 4.0 and Kp = 0.55 for 



n = 0.4, and a = 5.1 and K„ 



0.76 for n = 0.8. 



0.76 for n — 0.4 and n — 0.8, respectively. The values 
of Kp shown in the inset of Fig. for small values of J 
and L = 40 should be seen as upper Hmits. From the 
uncertainty of about 0.05 in the values of Kp for fixed 
L, we estimate the error in the extrapolated values to be 
<0.1. 

Now that we have accurately determined the exponent 
Kp we should be able to describe all correlation func- 
tions, since the only parameters needed are Kp and the 
Fermi momentum kp {K^ = 1 because of SU(2) symme- 
try). In particular, we can cross-check our results with 
the spin-spin correlation function (Eq.EJ. Previous work 
showed that for small values of the Kondo coupling J 
the size of the Fermi surface is smallfS^ so that kp is 
fixed (see also Ref. 0). To eliminate the effect of the 
open boundaries on the spin-spin correlation function, 
we considered a large system and averaged the correla- 
tions over pairs of sites separated by a given distance j to 
get (S-^ (0) • S"^ (/)), as discussed by other authorsiSii^iifi 
In Fig. 13 the circles correspond to this averaged spin- 
spin correlation function P{1) — (S^ (0) • S-^ i^)) aveg ^'-"^ 
n = 0.4 and n = 0.8 at J = 0.35 and L = 120. We 
restricted the values to the interval 30 < / < 90, because 
the TLL description only makes sense asymptotically and 
large values of I may be compromised by the open bound- 
aries. A direct attempt at fitting the data of Fig. with 



Is Kp < 0, which is clearly incorrect. We be- 



lieve the discrepancy is due to logarithmic corrections, 
which are well established in other models with SU(2) 
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Figure 6: Spin gap as function of 1/L for the densities n = 
0.5, n = 0.8 and n = 1 (J = 1.2). The data for n = 1 have 
been multiplied by 10~^ for comparison. 



symmetry, e. g. the Heisenberg modeli^fl^ Thus, as- 
suming a generic form 



(S^(O).S^(x)) 



1 



{ttxY 



Bi 10 



we can produce an excellent fit to the numerical results 
using the exponents Kp independently obtained before and 
only two fitting parameters: a = 4.0 {a = 5.1) and 
Bi = 0.17 {Bi = 0.11) for n = 0.4 {n = 0.8). This 
is seen as the dashed line in Fig. Logarithmic correc- 
tions thus appear to be crucial for a complete description 
of the spin correlations of the one-dimensional Kondo lat- 
tice model. As far as we know, this point has not been 
stressed before. We note that the above values of a dif- 
fer from the expected values of l/2i^ A more accurate 
determination of the exponent of the log correction may 
require much larger system sizes. 

Finally, we would like to address the quarter-filled 
case, n = 1/2, which has been shown to exhibit spin 
dimerization.? At this fiUing, the system has a charge gap 
and the charge sector cannot be described as a TLL.- Fur- 
thermore, since the spins are dimerized, we would naively 
expect a finite spin gap, as in the frustrated Ji — J2 
Heisenberg modelf22i We would now like to show that in 
fact the spin sector is gapless. In Fig. El we show the spin 
gap as a function of the lattice size L for J = 1.2 and den- 
sities n — 0.5, n — 0.8, and n — 1.0. We chose to work 
with a large J value in order to produce large spin gaps, 
since the spin gap generally increases with J. At half- 
filling, where the system is known to be fully gappedji^ 
the data clearly tend to saturate at a non-zero value in 
the thermodynamic Hmit. By contrast, at n = 0.5 and 
n = 0.8 the data strongly indicate that the spin sector is 
gapless. Thus, the dimerized phase discovered in Ref. 
has a charge gap but no spin gap. In that reference, the 
effect of dimerization in the localized spin sub-system on 
the conduction electrons was discussed. If we integrate 



5 



out the local moments, an effective exchange interaction 
among the conduction electrons is generated. This is in a 
sense the "complement" of the RKKY interaction, which 
induces an effective exchange interaction between local 
moments once the conduction electrons are integrated 
out. This effective exchange interaction is proportional 
to the static spin susceptibility of the dimerized locaHzed 
spins 

^fe//~ J'^x'a-A:)s, -Sfc. (11) 

jk 

If only nearest neighbor terms are retained — k) = 
5j^k+iD{j), where D{j) {—lyDo is the dimer order 
parameter. This leads to a staggered exchange interaction 
between conduction electron spin densities 

Heff-Y.(-iy^r^,+i- (12) 

This kind of interaction can be analyzed through 
bosonizationpifi Among the many terms that are 
generated, those which involve combinations like 
e-i'^'^!'^^p\^{x)iPLix)'iljli{x + l)iJL{x + l) will have just the 
right oscillating factor to cancel the (—1)-' in Eq. 
since Akp = tt at quarter filling. One of the terms is 



^ sin(2-v/2(/'p) (in the notation of R.ef . . which is rel- 
evant and opens a charge gap if i^p < 1^ The DMRG 
results show that this condition is fulfilled throughout 
the phase diagram. Thus, the above analysis seems to 
be consistent with the presence of a charge gap. A term 
of the form ~ sin{2^/2(f> cr) is also generated. However, it 
is marginal if K^- = 1 and only generates a spin gap if 
the coefficient has the right sign. We conclude that the 
presence of a charge gap and the absence of a spin gap 
we find is consistent with the bosonization analysis. 

In conclusion, we have presented a systematic study 
of the non-universal exponent Kp in the Kondo lattice 
model, as a function of the conduction electron density 
and the Kondo coupling. The qualitative behavior of the 
charge structure factor C (g) in the weak and strong cou- 
pling limits could be ascribed to free spin-i and spinless 
fermions, respectively. We also showed that the spin cor- 
relation function can be described within a Tomonaga- 
Luttinger liquid scheme only if logarithmic corrections 
are included. Finally, we have demonstrated that, al- 
though the charge sector has a gap at quarter fiUing, 
there are gapless spin excitations. 
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